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I .   INTRODUCTION 

Several  experimental  situations  involve  administering  a  stimulus 
to  an  experimental  unit  and  then  observing  a  binary  response.  Consider 
the  following  case  :  An  experimenter  can  apply  a  stimulus  at  various 
levels.  Subjects  are  assigned  at  random  to  different  levels  of  the 
stimulus  and  for  each  subject,  a  binary  response  is  recorded.  From 
this  set  of  observations,  one  determines  a  model  to  describe  the 
relationship  between  the  stimulus  and  the  proportion  responding  to  the 
stimulus.  The  model  is  called  a  dose-response  model.  These  types  of 
models  can  be  applied  to  a  lot  of  practical  problems  such  as  those 
presented  by  Chaterjee(1977) ,  Finney(1971) ,  and  Milliken(1982) .  They 
are  notedly  common  in  the  field  of  pharmacology,  in  bioassay,  where 
the  levels  of  the  stimulus  may  represent  different  doses  of  a  drug  or 
pesticide  and  the  binary  response  is  either  death  or  survival. 

Suppose  a  pesticide  is  tested  at  k  different  levels.  At  the 
level  of  dosage  X. ,  let  m.  be  the  number  of  insects  that  die  out  of  a 

total  n.  exposed  to  the  pesticide  (i  =-  1,  2,...,k).  We  want  to 

estimate  the  relationship  between  dose,  X.,  and  the  proportion 

dying,  P. . 

The  stimulus-response  relationships  have  generally  been  found  to 
be  nonlinear.  One  nonlinear  function  which  has  been  used  to 
approximate  the  relationship  between  the  dose,  X.  and  the  response 

probability,  P   is 


+  fix*L 


-       o  ^1  1 
1  +  e 


(  1.1  ) 


where  0     and  0X      are  unknown  parameters. 

This  relationship  is  referred  to  as  the  logistic  response 
function  and  has  the  shape  shown  in  Fig.  1.1.  The  corresponding 
analysis  is  referred  to  in  the  literature  as  logit  analysis.  The 
logistic  function  is  bounded  between  0  and  1,  and  is  monotonic. 

Two  estimators  of  0     and  0.    in  the  model  have  been  proposed  in 

the  literature.  One  estimator  is  the  maximum  likelihood  estimator 
(MLE)  and  the  other  is  the  method  proposed  by  Berkson  in  1953.  This 
study  aims  to  further  investigate  and  compare  the  behaviors  of  the  two 
estimators  with  respect  to  some  dose-response  designs  by  computer 
simulation. 
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II.   DERIVATIONS 

As  earlier  noted,  the  logistic  response  function 

fi   +  a.x. 

Pi"  ■  °    l   l  , 


1  +  e 


is  commonly  used  Co  approximate  the  relationship  between  the  dose  X. 

and  the  response  probability  P  .  The  main  concern  is  to  estimate  the 

parameters  j0   and  0. .  Two  estimators  of  the  logistic  regression 

coefficients  have  been  proposed  in  the  literature  and  are  derived  in 
this  section. 

A.  Maximum  Likelihood  Estimator 

Let  X.-  ith  dose  i-  1,2 k  levels, 

n.—  number  of  experimental  units  exposed  to  the  stimulus, 
m.-  number  of  experimental  units  reacting  to  the  stimulus, 

p.-  _i  -  observed  'mortality  rate' 

ni 

P.-  'true  mortality  rate'  and 

Q.-  1  -  P.. 

l       l 


It  is  usually  reasonable  to  assume  that  m.  is  a  Binomialfn. ,  P.) 

l  11 

random  variable.  Thus,  the  probability  density  function  of  m.  is 


for  I  -  1,. ,k. 

The  log  likelihood  function  of  the  model  is 


k    r   .    k  k 

In  L  -2  In   ni   +.2.m, In  P.   +  .2. (n. -m, )ln  Q.  .   (2.2 
i-1    I  m  J    i-l  l    i    i-l   l  i     l 


where  P.  -  e ,  Q.-  1-P.  . 

1        8  +  B.X.  1 

1  +  e  °   X  x 

Differentiating  In  L  with  respect  to  8     and  8,    and  setting  the 

corresponding  derivatives  equal  to  zero,  we  obtain  the  normal 
equations 


k 
i-l  V  pi  "  Pi  )   ~  °  (  2.3  ) 

and 

k 
l5l  niV  Pi  "  Pi  )   '  °  (  2.4  ) 


Solving  the  normal  equations  by  iteration  gives  the  ML  estimators  for 

8      and  fln . 
'o     ^1 


B.  Berkson's  Method 


Let  X..  m. .  n.  ,  p.,  q.,  P.,  Q.  denote  the  same  quantities  as 
1   1   i  ri  ni   1   l  * 

in  Section  II,  A.  In  the  following,  we  outline  the  derivation  of 


Berkson's  estimators  of  8     and 


It  is  well-known  that 


AN 


P.  , 
l 


Pi(1-V 


(2.5) 


where  "-AN"  means  "asymptotically  normally  distributed" .  Then 


In 


I  l-PtJ 


-  AN 


+  fl,X 


1  1"  n.P.(l-P.) 
l  l    l 


(2.6) 


(  See  for  example  p. 118,  Theorem  A,  Serf ling  (1980)  ).  Let 

for  1-1,..,  k. 


1.  -  In 


pi 


Then  one  can  write: 


1,  -  8     +  0.X,  +  e, 


where  e.  -  N(0,  (n.P.(l-P.))   )     for  1-1,2, ...k, 
1  i  1    l  lit  i 


and  e  x  ,  .  .  ,  e,   are  independent. 


(2.7) 


One  first  obtains  the  weighted  least  squares  estimators  of  0      and 
3X    by  using  (2.7)  under  the  assumption  that  P.'s  are  known.  Then 


replacing  P.  by  p,  in  the  weighted  least  squares  estimators  of  0     and 
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fl-      one  obtains  Berkson's   estimators   of  6     and  B^    as   follows: 
^1,  'o  ^1 

k  ,    „2  .2,    n.p.q. 1.  .2.    n.p.q.X. 

.2,    n.p.q.l.X.    -        l-l      i,ri    l    l  l-l      i*i^i    l 

l-l      1*1^1   i   l  k  '       

a  .  2.,    n .  p  .  q .  (2.8) 

„  l-l      ir  li 

"l      ~  k  2 

k  _  2  (.2,    n.p.q.X.) 

.2,    n.p.q.X.  -  vi-l.     i*iMi   i/ 

l-l      iri^i   l  fc-^6 

.2,    n.p.q. 
i-1      iri^i 


and 


^    -  Ji  "i?!"!1!   -    3i  Ji  yt^i  <■  2-9  > 


.2.    n.p.q. 
1-1      lrlni 


III.   REVIEW  OF  LITERATURE 

Other  than  the  logistic  model  discussed  in  Chapter  1,  an 
alternative  model  in  which  the  response  function  is  represented  by  the 
cumulative  normal  distribution  function  has  also  been  proposed.  The 
cumulative  distribution  function  of  the  normal  distribution  has  a 
shape  similar  to  that  of  the  logistic  function  (see  Fig  1.1).  The 
resulting  model  is  called  the  probit  model  and  is  extensively  covered 
by  Finney(1971) . 

Milliken(1982)  discussed  dose-response  analyses  in  the  light  of 
tolerance  distributions.  Tolerance  of  the  experimental  unit  defined  as 
the  level  of  the  stimulus  at  which  a  response  first  occurs,  he 
described  the  tolerance  distribution  as  the  collection  of  tolerances 
for  a  population  of  experimental  units  and  that  the  assumed  form  of 
the  tolerance  distribution  determines  the  type  of  analysis.  Aside  from 
the  most  commonly  used  probit  and  logit  analysis  which  uses  the  normal 
and  logistic  curves  as  tolerance  distributions,  respectively,  models 
resulting  from  the  use  of  exponential  and  e  to  the  e  function  as 
tolerance  distributions  were  also  presented. 

For  the  logit  analysis,  several  studies  have  been  undertaken 
since  the  maximum  likelihood  estimator  (MLE)  and  Berkson's  method 
appeared  in  the  literature.  The  MLE  was  first  advocated  for  general 
application  by  R.A.  Fisher.  For  this  particular  model  though,  the 
solution  requires  an  iterative  procedure  as  opposed  to  a 
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computationally  simple  estimator  suggested,  by  Berkson(1953) .  This  is 
in  fact  one  of  the  factors  which  motivated  this  study  of  comparison. 

Further  interest  was  brought  about  by  some  notable  results 
presented,  in  a  theoretical  study  by  Davis(1985).  One  of  the 
interesting  implications  of  her  results  was  that  Berkson's  method 
should  not  be  used  when  the  number  of  design  points  is  large  and  the 
average  number  of  observations  per  design  point  is  small.  Furthermore, 
she  pointed  out  that  when  the  number  of  design  points  is  large, 
Berkson's  estimator  is  consistent  if  the  average  number  of 
observations  per  design  point  is  large.  These  results,  which  will 
eventually  (in  Chapter  V)  be  compared  with  the  outcome  of  the 
simulation  study,  further  motivated  this  area  of  research. 


For  the  logit  transformation,  the  value  of  In 


defined  for  p.-  0  or  1.  Chatterjee  and  Price(1977)  noted  that  for 

fitting  a  dose-response  relationship,  such  points  are  usually  omitted, 
as  there  is  considerable  degree  of  uncertainty  about  the  exact  dosage 
at  these  response  levels.  Berkson's  proposal(1955) ,  which  was  used  in 
this  study,  was  to  use  the  following  rule: 


P,- 


2n. 
l 

1 


m.  -  0 


1 
2n. 


where  p.,  n. ,  and  m.  are  as  defined  in  Chapter  II. 
ri   li  K 


IV .   METHODOLOGY 

The  special  case  of  a  simple  logit  model  with  one  independent 
variable  was  examined.  For  all  experiments,  the  parameters  of  the 
logistic  function  were  0  -   -4.12  and  fi  -   1.18.  Nine  dose-response 

designs  were  investigated:  three  five-dose  designs,  three  ten-dose 
designs  and  three  fifteen-dose  designs.  The  experiments  were  conducted 
with  an  equal  number  of  observations  per  cell,  n  -n  for  all  i.  In  all 

three  main  simulations,  the  cell  sizes  were  n  -  10,  20,  30.  Therefore, 
the  total  sample  sizes  for  the  five-dose  designs  were  N  -  50,  100, 
150;  N  -  100,  200,  300  for  the  ten-dose  designs  and  N  -  150 ,  300,  450 
for  the  fifteen-dose  designs. 

For  each  of  the  9  simulations,  500  sets  of  data  were  generated 
through  GAUSS(1987)  (  for  the  algorithm,  see  Appendix  1)  and  were 
analyzed  using  the  two  estimation  methods.  Here,  for  the 
aforementioned  P's(fi   — 4.12  and  0.-1.18)  and  specified  X.'s(k  doses 

V  "lXi 
equally  spaced  on  log  doses  1-6),  P.  -  e were  computed 

1  0  +  A,X. 
i  .  O  'li 
1  +  e 

and  m.'s  were  generated  from  a  Binomial (n. ,  P.)  random  variable  where 

n. 's  were  the  specified  number  of  observations  per  dose  in  the  design. 

The  Berkson  estimates  were  obtained  using  a  GAUSS  algorithm  while  that 
of  maximum  likelihood  were  computed  using  PROC  NLIN  of  SAS(1982). 
As  a  basis  for  comparison,  biases,  variances  and  mean  square 
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errors(MSE)  of  the  estimated  regression  coefficients  (  0    ,    0.       )  were 
calculated  for  each  design.  In  addition,  as  a  measure  of  how  the 
estimated  p's(  p  )  compare  with  the  actual  P's,  the  quantity 


A  (  pt  -  pi  > 


was  computed  for  each  dose-response  data  set  and  averaged  over  500 
data  sets  under  each  design.  This  in  essence,  measures  on  the  average, 
how  close  the  predicted  logistic  regression  curve  is  from  the  actual. 
Throughout  the  report,  the  quantity 


j-si 


X  <p-  -  p->2 

l-l  *i    i 


(  4.1  ) 


where  s  -  number  of  data  sets(in  this  case  500) 
will  be  referred  to  as  MSE(p) . 
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V.   RESULTS  AND  DISCUSSION 

A.  Maximum  Likelihood  vs.  Berkson's 

Simulation  estimates  of  the  bias,  variance  and  MSE  of  the 

sampling  distributions  of  the  estimated  B     and  B.    are  presented  in 

Table  5.1.  MSE(p)(see  4.1)  are  also  included  in  the  table. 

1 .  Biases  of  3     and  3, 

lo n. 

Simulation  results  indicate  that  the  Berkson  statistic  is  more 
biased  than  the  maximum  likelihood  estimator  (MLE) .  For  the  small 
sample  size  within  dose  designs  (k  -  5 ,  10,  15;  n-10)  the  differences 
were  not  great.  Berkson  biases  were  only  .78  to  2.5  times  larger  than 
MLE's.  These  are  graphically  demonstrated  in  Fig.  5.1-5.12. 
Differences  were  most  pronounced  in  larger  sample  size  designs  (n-20, 
30  under  all  k) ,  as  illustrated  in  Fig.  5.13-5.32.  In  these  cases,  MLE 
essentially  dominated  Berkson's  (with  Berkson  bias  being  as  much  as 
176  times  that  of  MLE)  except  in  the  5  dose,  n-30  design.  Here,  MLE 
bias  was  about  1.5  times  larger  than  Berkson's.  This  was  the  only  case 
where  Berkson's  fared  better  than  the  MLE  in  terms  of  bias.  Figures 
5.33-5.36  suggest  though  that  Berkson's  advantage  is  small. 

2.  Mean  Square  Errors  of  8     and  3- 

Mean  square  error  results  for  the  regression  coeficients(/9  ,  0.  ) 

were  consistent  with  the  performance  of  the  estimators  with  respect  to 
biases.  In  this  case  though,  the  discrepancies  were  rather  small. 
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Berkson  MSE's  were  only  7%  to  56%  higher  than  MLE's,  except  again  for 
the  5  dose,  n-30  design  where  MSE  of  maximum  likelihood  estimators  for 
8     and  B^   were  22%  and  38%  more  than  their  Berkson  counterparts. 

3.  Variances  of  8     and  fl, 
_o ^1 

The  variances,  on  the  other  hand,  exhibited  a  slightly  different 

trend.  Though  not  extremely  nor  uniformly  better,  the  Berkson 

estimates  showed  less  variability  as  seen  in  Table  5.1.  This  was 

exhibited  in  all  of  the  5  dose  designs  (n-10,  20,  30)  and  in  the  15 

dose,  n-20  design  where  MLE  variances  were  1  to  30%  larger  than  the 

corresponding  Berkson' s.  Plots  in  Fig.  5.29-5.32  particularly  reveal 

this  Berkson's  rather  slight  edge  over  MLE  in  some  designs.  However, 

maximum  likelihood  still  yielded  slightly  smaller  variances  for  both 

0o   and  01   in  the  15  dose,  n-30  design(Fig.  5.13-5.16).  In  all  other 

cases  the  trend  was  inconclusive.  There  was  inconsistency  in  the 
behavior  of  the  variances  within  a  design.  That  is,  either  Berkson's 

gave  lower  variance  for  8     and  higher  variance  for  8     or  vice-versa. 

4.  Mean  Square  Error  of  p 

In  most  cases,  experimenters  are  not  particularly  interested  in 
the  specific  values  of  8     and  8     but  are  concerned  with  estimating  P., 

which  is  the  actual  proportion  of  experimental  units  reacting  to  the 
stimulus  at  a  certain  level  of  dosage.  A  quantity  which  measures  the 

'combined  ability'  of  8     and  8    ,  in  predicting  p.  was  computed.  This 
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was  the  MSE(p)  discussed  in  Chap.  IV(4.1).  It  essentially  measures  how 
close,  on  the  average,  the  predicted  logistic  regression  curve  is  to 
the  actual  curve . 

The  maximum  likelihood  estimator  showed  complete  superiority  in 
this  respect  as  evidenced  in  Fig.  5.37-5.54.  Berkson's  estimates 

produced  MSE(p)'s  which  were  3%  to  69%  higher  than  the  MLE 
counterparts.  It  was  also  observed  that  the  discrepancy  increased  as 
the  sample  size  increased  within  a  specified  number  of  doses  (k) . 

B.  Berkson  Estimator  in  9  designs 

For  a  fixed  number  of  doses (k) ,  as  the  sample  size  increased  the 
bias  decreased,  the  variances  were  reduced,  and  MSE's  decreased  in  all 
of  the  designs  as  anticipated.  It  is  again  worth  noting  that  the 
Berkson  estimator  was  at  its  best  as  far  as  bias  was  concerned  in  the 
5  dose,  30  observations  per  dose  design. 

C.  Maximum  Likelihood  Estimator  in  9  designs 

As  the  sample  size  n  increased  within  a  fixed  number  of  doses(k), 
the  improvement  in  bias,  variance,  and  MSE,  that  is,  in  the 
performance  of  the  estimator  in  general,  was  not  as  great  as  that  in 
Berkson.  This  shows  the  relative  stability  of  the  MLE  regardless  of 
the  design.  In  fact,  for  k-5  and  k-10 ,  biases  were  even  smaller  when 
n-20  than  when  n-30. 

D.  Comparison  of  Designs  with  Equal  Total  Sample  Sizes  N 

This  part  of  the  analysis  is  focused  on  the  practical  aspect  of 
experimentation.  Most  often,  experimenters  are  faced  with  a  fixed 
number  of  experimental  units  at  hand  and  the  immediate  problem  is 
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choosing  a  proper  dose-response  design  to  most  closely  predict  p 
based  on  a  fitted  logistic  model.  Is  one  better  off  using  more  design 
points  or  using  more  observations  per  design  point? 

1.  N-100  :  k-5  n-20  vs.   k-10  n-10 

The  5  dose,  20  observations  per  dose  design  resulted  in  lower 

MSE(p)  for  both  methods  as  shown  in  Table  5.1  and  Fig.  5.39-5.40, 
5.43-5.44.  No  definite  trend  was  discerned  among  the  0    's  and  f)   '  s  not 

only  in  this  case  but  in  the  other  two  comparisons  below  as  well. 

2.  N-150  :  k-5  n-30  vs.   k-15  n-10 

Here,  Table  5.1  and  Fig.  5.41-5.42,  5.49-5.50  conspicuously 
reveal  the  superiority  of  the  5  dose,  30  observations  per  dose  design 
over  the  15  dose,  10  observations  per  dose  design,  in  predicting  the 
proportion(p)  of  experimental  units  responding  to  the  stimulus. 

3.  N-300     k-10  n-30  vs.   k-15  n-20 

The  15-dose,20  observations  per  dose  design  is  slightly  better  in 

predicting  p  as  evidenced  by  its  low  MSE(p)  in  both  designs.  Fig. 
5.47-5.48,  5.51-5.52  attest  to  this. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

These  simulation  results  concur  with  the  conclusion  by 
Davis(1985)  that  Berkson's  method  is  not  advisable  when  the  number  of 
design  points (k)  is  large  but  the  average  number  of  observations (n) 
per  design  point  is  small. 

The  maximum  likelihood  demonstrated  superiority  in  every 
criterion  except  for  the  variances  under  all  designs  but  one (5  dose, 
30  observations  per  dose  case).  Most  importantly,  it  fared 

consistently  better  in  the  MSE(p)  criterion  which  is  probably  of 
utmost  consideration  in  practical  applications. 

Although  the  maximum  likelihood  is  the  method  recommended  for  use 
based  on  the  results  of  this  study,  the  simulation  results  did  show 
some  indication  that  Berkson's  method  may  be  better  for  experiments 
when  the  number  of  design  points  is  small  but  the  number  of 
observations  per  design  point  is  large.  Further  studies  may  be  in 
order  regarding  that  situation. 
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VIII .   APPENDIX 

A. Sample  Gauss  Program  :  Generation  Of  Data  Sets 

seed-45989; 
1-451; 

do  while  1  <-500; 
dose-1; 

do  while  dose  <-  6 ; 

pow— 4.12  +  (1.18*dose); 
pt-l/(l+exp(-pow)) ; 
r annum— rndus ( 20 , 1 , seed) ; 
dead-r annum [ . , 1] .<  pt; 

dnum-sumc(dead) ;  if  dnum— 0;dnum-20*(l/40)  ; 
endif; 

if  dnum—  20;dnum-20*(l- (1/40) )  ; 
endif; 

if  dose  —  1  ;  deaths— dnum;x-dose;  totl-20 ; 
ptrue— pt   ;  sampleno=i ; 
else;  deaths  - 

deaths | dnum ; x-x | dose ; totl-totl | 20 ; 
ptrue-ptrue |pt ; sampleno-sampleno | i ; 
endif; 
dose-dose+ . 555 ; 
endo  ; 


/*  to  form  the  dose  response  matrix  x  */ 


dataset  -  x~totl-deaths~ptrue-sampleno; 

format  /ml/rd  11,8; 

output  file-datl020j  on;  dataset; 

i-i+1 ; seed-seed+100 ; 

endo; 
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B.  Sample  Gauss  Program  :  Berkson  Estimates 

/*  Berkson  :  this  program  computes  Berkson  estimates  */ 

p-deaths . /totl ; 

pprime-ln(p./(l-p)); 

weight-totl.*p.*(l-p) ; 

wpprimex-weight . *pprime . *x ; 

wprim-weight . *pprime ; 

wx-weight.*x; 

wxsqr-weight.*(x.*x)  ; 

bl-( (sumc(wpprimex) ) - ( ( (sumc(wprim) )*(sumc(wx) ) >/(sumc (weight) ) ) )/ 

((sumc(wxsqr))-(((sumc(wx))A2)/(sumc(weight)))) ; 
bo-(meanc(pprime))-(bl*(meanc(x))) ; 

/*  to  compute  the  mse  */ 

phat-(exp(bo+(bl*x)))./d+(exp(bo+(bl*x)))); 
ssep-meanc( (phat-ptrue) A2) ; 
ssebo-(bo+4.12)A2; 
ssebl-(bl-1.18)A2; 

/*  to  output  into  a  data  set  */ 

if  i— 1  ;  matbl-bl;matbo— bo;matp-ssep;matsbl-ssebl;matsbo-ssebo; 
else; 
matbl-matbl|bl;matbo-matbo|bo;matp-matp | ssep;matsbl-matsbl | ssebl; 
matsbo-raatsbo | ssebo; 
endif ; 
seed-seed+100 ; 
i-i+1; 

endo; 

/*  to  print  the  matrix  */ 

e  s  t-matb l~ma tbo-matp~matsbl~matsbo ; 

format  /ml/rd  11,8; 

exvalue-meanc(est) jdev-stdc(est) ; summary-exvalue-dev ; 

lprint  "ESTIMATES  AND  MEAN  SQUARE  ERRORS(CHIT1020) "  est; 

output  file-b;chitl020  on; 

lprint  "SUMMARY  STATS(chitl020) "  summary'; 
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C.  SAS  Program  :  Maximum  Likelihood  Estimates 


*  THIS  PROGRAM  COMPUTES  THE  * 

*  MAXIMUM  LIKELIHOOD  ESTIMATES  * 

*  FOR  THE  LOGIT  MODEL  * 

* - *; 

%MACRO  ANALYSIS; 
%DO  1-1  %T0  500; 

DATA  ANA&I;SET  RJR.SIM1520; 

IF  SAMPLENO-&I; 
DATA  ANA26J;SET  ANA&I ; 

COUNT-DEATHS ;  Y-l ;  OUTPUT ; 
COUNT-TOTAL-DEATHS  ;  Y-0  ;  OUTPUT ; 
PROC  NLIN  NOHALVE  SIGSQ-1   /*ALGORITHM  FROM  SAS82  STATISTICS 

P36*/; 

PARMS  BO— 4.6  Bl-1 . 3/* INITIAL  ESTIMATES  OF  BO  AND  Bl  FROM 
PROC  REG*/; 

E-EXP(BO  +  B1*D0SE); 

P-E/(l+E); 

MODEL  Y-P; 

W-1/(P*(1-P)); 

_WEIGHT_-W*COUNT ; 

DER-E/((1+E)**2); 

DER.BO-DER; 

DER.B1-DER*D0SE; 

OUTPUT  OUT-RESULTS  P-PHAT  PARMS-BO  Bl ; 
DATA  PREFINAL; 

KEEP  BO  Bl  SSE; 

SET  RESULTS; 

SSE-((PHAT-PTRUE)**2) ; 
PROC  MEANS  MEAN  NOPRINT; 

OUTPUT  OUT-FINAL  MEAN-BO  Bl  MSE; 
PROC  APPEND  BASE-FINALDAT  DATA-FINAL; 
%END; 
%MEND; 
%ANALYSIS; 
DATA  RJR1.ML1520;SET  FINALDAT; 
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ABSTRACT 

Two  estimators (Maximum  Likelihood  and  Berkson' s(1953))  of  the 
logistic  regression  model  have  been  proposed  in  the  literature.  The 
two  estimators  were  compared  by  simulation  under  9  different  dose- 
response  designs.  Maximum  likelihood  emerged  as  the  better  estimator 
in  all  cases  considered  except  one.  Berkson' s  method  showed  some 
promise  in  a  few  design  point-many  observations  per  design  point  set- 
up. 


